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Abstract 

Coil-stretch hysteresis in extensional flows of dilute polymer solutions is shown to be strongly en- 
hanced as concentration c increases towards c* , the critical overlap concentration. Large transverse 
fluctuations within the window of extension rates where hysteresis occurs cause greater intermolec- 
ular overlap between partially stretched molecules than at equilibrium. Hydrodynamic screening 
thus sets in at values of c/c* « 1 for such chains, increasing their mean frictional coefficient and 
allowing for stretched states to be sustained at very low extension rates. Hysteresis progressively 
diminishes beyond c* , and vanishes altogether when the crossover from Zimm to Rouse dynamics 
is complete. The concentration dependence of the width of the hysteresis window is thus non- 
monotonic, with a large maximum at c/c* = 1. 



The dynamics of flexible polymers in a dilute solution are determined by the interplay 
between forces arising from thermal fluctuations, and intramolecular free-energy and hy- 
drodynamic interactions. This is embodied in the relationship Az = Cz/4i?o between the 
characteristic "Zimm" relaxation time Az, the friction coefficient £z of isolated polymer 
coils at equilibrium and the average stiffness H of the internal resistance of equilibrium 
coils to stretching. The latter in turn is related to the equilibrium end-to-end distance R : 
H = 3k B T/Rl, where T is the absolute temperature of the solution and k B denotes the Boltz- 
mann constant. The seminal work of De Gennes [1], Hinch [2] and Tanner [3] suggested that 
in flows with a dominant extensional component, a sharp coil-to-stretch transition occurs for 
long molecules when the extensional rate e is such that the Weissenberg number Wi = e\z 
exceeds a critical value Wi c _ s = 1/2, and molecules unravel and stretch significantly in the 
flow. The existence of this transition has been beautifully demonstrated in single-molecule 
experiments [4] and is now thought to play an important role in a number of complex flows 
of dilute polymer solutions [5; 6]. 

Both the molecular stiffness H and friction coefficient ( are in general functions of molec- 
ular conformation. In a solution close to its 0-state where the effects of excluded volume 
interactions are negligible, molecular resistance to stretching has entropic origins. Although 
for small deformations this entropic resistance is linear in molecular stretch and the stiffness 
H ~ H , it diverges as the end-to-end distance R approaches the contour length L [7; 8]. 
Intramolecular hydrodynamic interactions result in an increase in ( with R, from at equi- 
librium towards the value given by Batchelor's result for isolated slender rods [9]. As a 
consequence of the non-linearities in conformation dependence of H and (, the coil-stretch 
transition in extensional flows is associated with strong hysteretic effects [1-3]. The steady 
state in an extensional flow is primarily the result of a balance between internal entropic 
resistance to stretching and frictional drag forces exerted by the flowing solvent, which ef- 
fectively occurs when (e = 2H. In very weak flows, this balance is achieved at R > R 0: and 
when Wi > Wi c _ s , only highly stretched states (R < L) are stable. Within a range of extension 
rates such that Wi s _ c < Wi < Wi c _ s however, both the coiled and stretched states are stable, 
and it is possible to kinetically trap molecules in either of these stable states depending on 
initial conditions [1]. The lower bound of the hysteresis window Wi s _ c can be interpreted 
as a critical strain rate for an ensemble of stretched molecules to transition to the coiled 
state. Single molecule experiments and simulations [10; 11] and rheological measurements 
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[12] have confirmed the existence of this window. 

It is known that in a dilute solution, the size of the hysteresis window measured in terms 
of the ratio Wi c _ s /Wi s _ c is proportional to ratio of £ s , the drag coefficient of a nearly fully 
stretched molecule, to the coiled-state friction (z- In a concentrated polymer solution on the 
other hand, it is expected that conformational dependence of the friction coefficient vanishes 
since interpenetration of molecules completely screens out solvent-mediated hydrodynamic 
interactions. The behaviour of any single molecule is "Rouse-like" as opposed to the "Zimm- 
like" dynamics of isolated chains in dilute solutions [8; 13]. Both coiled and stretched 
molecules have the same Rouse friction coefficient Cr in a concentrated solution, and hence 
no coil-stretch hysteresis is expected in this case. 

The objective of this Letter is to show that as the concentration c is increased past the 
critical value c* at which equilibrium coils begin to overlap, the change in the width of 
the coil-stretch hysteresis window with the ratio c/c* may not be monotonic, but attains a 
maximum when c/c* ~ 1. The strong enhancement of hysteresis with concentration occurs 
for c/c* « 1, that is, at concentrations that would conventionally place solutions well in the 
dilute regime. 

Recent experiments [16] and molecular simulations have demonstrated that the behaviour 
of even very dilute polymer solutions in strongly stretching flows can exhibit non-trivial con- 
centration dependence. Some insight into the origin of this "self-concentration" phenomenon 
can be obtained from single-molecule simulations. Conventionally, the critical overlap con- 
centration i defined as c* = d^ 3 , where do is the size of isotropic coils at equilibrium, estimated 
here as d = i? /v3- Interpreting c as a number density of polymer molecules in solution, 
the fractional volume pervaded by isotropic coils at equilibrium is cd^ = c/c*. In extensional 
flows within the coil-stretch hysteresis window, conformational fluctuations of size d trans- 
verse to the principal stretching axis are large and comparable in magnitude to equilibrium 
fluctuations d since Wi < 1 (Fig. 1 a) within the window. Therefore, the pervaded volume 
id 2 of a molecule in its stable, stretched state with an average length £ » R along the 
principal direction of extension is much larger than at equilibrium (Fig. 1 b). The volume 
fraction pervaded by stretched molecules, 

Thus, when aspect ratio £/d » 1 but with transverse fluctuations d/d ~ 0(1), a solution may 
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FIG. 1. (a) Results of Brownian-dynamics simulations [14; 15] for gyration tensor component 
G rr of isolated chains transverse to the direction of stretching in steady uniaxial extensional flows; 
(b) mean pervaded volume in simulations of isolated linear bead-spring chains: the volume V is 
estimated from components of the gyration tensor G as V = G rr \/G zz , and Vq is the equilibrium 
volume of isotropic coils. Coil-stretch hysteresis is observed for Wi < 0.5; chains initially at equi- 
librium (A) and chains initially stretched to 90% maximum permissible length (▼) have different 
quasi-steady-states. The horizontal dashed line in (b) indicates the volume estimate assuming 
G zz = L 2 /12 (valid for rods of length L) and G rr = G rr ^/2. 



initially be very dilute at equilibrium with c/c* « 1, but <fi could be significantly larger in an 
extensional flow in the vicinity of the coil-stretch transition. Intermolecular interactions can 
be more significant than at equilibrium and it is even possible that partially stretched coils 
may overlap, as observed recently in multi-chain simulations [17]. The question is: how does 
self-concentration affect the dependence of the friction coefficient C on molecular stretch? 

As already mentioned, in the limit of infinite dilution, the rod-like friction coefficient 
Cr of stretched molecules of length £ and diameter d can be estimated from the theory for 
suspensions of slender rods [9] as Cr/Cz ~ (( / d ) / \n(£ / d) . For non-zero c, the average distance 
between stretched molecules h = (c£) -1 / 2 , and when d/h = 1 / 2 « 1, far-field hydrodynamic 
interactions modify the friction coefficient, and approximately for non-overlapping chains 
when the instantaneous <fi « 1 [9], 



Cr 



1 



(2) 



Cz \n(W-Mi + W d o' 

To the best of this author's knowledge, a theory for a solution of anisotropic chains near and 
beyond overlap has not been developed. We can nevertheless draw upon current understand- 
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FIG. 2. Schematic of length-scales in a solution of partially unravelled, overlapping polymer chains 
aligned along the principal stretching direction in a uniaxial extensional flow: each chain consists 
of (a) Nd = i\d domains of size determined by transverse fluctuations, and (b) correlation blobs 
of size ^h- 

ing of the Zimm-to- Rouse crossover for isotropic coils [8; 13; 18]. In a solution consisting of 
partially stretched chains all oriented along the axis of extension, progressive screening of 
hydrodynamic interactions with increasing overlap is described in terms of a hydrodynamic 
screening length £ h < d such that intramolecular hydrodynamic interactions persist only 
at length scales smaller than £ h , but are screened at larger length scales. Therefore, each 
molecule on average acts as a Rouse chain of N h "blobs", each of size £h (Fig. 2) [13; 18]. The 
quantities £h an d are determined by firstly noting that £h is the scale at which segments 
in any chain first encounter those from neighbouring chains; hence c£jjiVh = 1. Secondly, 
assuming that chain segments at length scales smaller than d retain a spatial distribution 
similar to the ideal random-walk at equilibrium, and £l/d 2 = (N K /N h )(N K / N d ) = N d /N h 
where N K = L 2 /Rl is the number of Kuhn segments in the molecule, and N d = i/d is the 
instantaneous aspect ratio (Fig. 2). Therefore, when > 1, £ h /d = <fi~ l and TVh = 4> 2 (i/d), 
and the total Rouse-like friction coefficient is 



In the regime of intermediate values of 0.01 < <fi < 1, a simple linear interpolation is used 
between the estimates given by Eqns. (2) and (3) above for Cr/Cz [15]. The friction of 
equilibrium coils ( c also increases with c/c* from the value for isolated coils Cz, firstly due to 




(3) 
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far- field hydrodynamic interactions when c/c* < 1, and then due to screening when c/c* > 1. 
For the sake of simplicity, changes Cc before overlap are neglected; that is, it is assumed that 
Cc = Cz when c/c* << 1. Beyond overlap, accounting for screening as outlined above gives 



Cz C* 



(4) 



Screening is considered to be complete when either ( T or Cc attain the value corresponding 
Cr. Therefore, they are limited in the model here to the maximum value set by Cr/Cz = 
(b K /d )N K = V3 N K , where b K = R\/L is the Kuhn length. Equation. (3) suggests that for 
stretched chains Cr = Cr at a concentration ct such that 



c 
c* 



-^(s)" (5)" ^ 

whereas for coiled molecules the Zimm-to- Rouse crossover is complete at a c^/c* = Cr/Cz = 
y/3N K . In other words, not only is the crossover for stretched chains complete at ctt « ct, 
but if I ~ L and (d/do) 2 ~ 0(1), the crossover occurs at c't/c* ~ iV K « 1, well within the 
dilute regime for solutions of long molecules. 
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FIG. 3. Variation of friction coefficent ( with molecular stretch I at fixed transverse chain 
dimensions (d/do) 2 = 0.5: the numbers alongside curves indicate c/c*. The maximum value for the 
friction coefficient Cr/Cz and the maximum permissible length £/do are both equal to \/3N K . 



The mean friction coefficient C is calculated here using a linear mixing rule to interpolate 
between Cr and Cc depending on polymer stretch [15]; Fig. 3 shows the variation of C with 
molecular stretch for c/c* values ranging from the limit of infinite dilution to nearly c^/c* 
when the transverse dimensions are held fixed such that (d/d ) 2 = 1/2. The friction coefficient 
Cs for hypothetical transversely bloated, fully-stretched chains with £ = L is much stronger 
than that of Cc when c < ctt. The ratio Cs/Cc becomes independent of concentration when 
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ctt < c < 1 since both ( s = Cr and ( c = Cz are constant. Above c*, the ratio Cs/d decreases 
with c as Cc "catches up" with ( s , and beyond ct, Cs/Cc = 1- 

Predictions for coil-stretch hysteresis are obtained with a mean-field conformation-tensor 
model whose dynamical equations for the evolution of the axial (zz) and transverse (rr) 
components of the conformation tensor M are [19]: 

dM zz AH Ak B T 

—j— =2eM zz -— M zz + — — , (6) 

dM rr . Ajf AH % , Ak B T 

= -eM rr -—M„ + ■ (7) 

The model consistently accounts for the effects of thermal fluctuations and flow, along 
with conformation-dependent molecular friction and entropic stiffness. Here, H = H Q (1 - 
l/iV K )[l - (M zz + 2M rr )/L 2 Y 1 is the nonlinear effective stiffness of the resistance to stretch- 
ing, and the friction coefficient C is calculated as described above [15] after identifying the 
instantaneous i and d of partially stretched coils with M X J Z and M*/ 2 , respectively. The near- 
equilibrium relaxation time Ao = ( c Rl/(12k B T), so that Ao = Az when c/c* < 1, and varies 
as Cc f° r concentrations above c*. The polymer contribution to the extensional viscosity, 
rj p = cH(M zz - M rr )/e. On rescaling model equations with R and A as the characteris- 
tic length and time scales, the sole dimensionless model parameters are N K , c/c* and the 
Weissenberg number, Wi = e\o. 

Unlike the results shown in Fig. 3, the hysteresis predictions in Fig. 4 are with the trans- 
verse size d/d of chains consistently determined at every Wi by the dynamics in Eq. (7). It is 
observed that for iV K = 5 x 10 3 , hysteresis windows widen first with increasing concentration 
in the dilute regime (Fig. 4 (a)) until about c/c* ~ 0.1, after which a decrease in window 
size commences only beyond c* , with hysteresis vanishing completely at c/c* ~ 100 (Fig. 4 
(b)). The non-monotonic concentration dependence of Wi c _ s /Wi s _ c is shown in Fig. 5, which 
compares well with the estimate Cs/Cc calculated a priori with the assumption of constant 
d/d = l/y/2. 

Although the discussion above is for solutions close to the 0-state, similar arguments 
are expected to hold when excluded-volume interactions are important. Indeed, the self- 
concentration effect may be marginally greater because the hysteresis window size at infinite 
dilution (horizontal line in Fig. 5) is smaller in good solvents and may even vanish completely 
[21; 22], but Cr (and hence the maximum in Fig. 5) is expected to be unaffected by solvent 
quality. An analysis of the dependence of the scaling exponents shown in Fig. 5 on solvent 




FIG. 4. Predictions of the conformation-tensor model for coil-stretch hysteresis in ry p , the polymer 
contribution to steady-state extensional viscosity (scaled by ck B TXq), shown for the sake of clarity 
separately for (a) c/c* < 1, and (b) c/c* > 1; the numbers alongside curves indicate c/c* values. 
Symbols shown in (a) are data for the window predicted by Brownian dynamics simulations (details 
shown separately in [15]). The thick gray curves in (a) are predictions obtained with the model 
in the infinite-dilution limit, while that in (b) is the prediction by the (FENE-P) model [20] with 
Rouse- like conformation- independent C = Cc- 
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FIG. 5. Concentration-dependence of the width of the coil-stretch hysteresis window: symbols 
are predictions of the conformation tensor model. The horizontal line is the width predicted in the 
infinite-dilution limit, while the curve bounding the symbols from above is estimated by the ratio 
Cs/Ccj w hh Cs being calculated from data in Fig. 3 for fully-stretched chains with large transverse 
fluctuations, that is £ = L, and (d/d^) 2 = 1/2. 

quality requires consideration of the changes in H with solvent quality [22; 23] which is 
beyond the scope of this Letter. 

There appears to be indirect experimental evidence for the predictions above from mea- 
surements of the characteristic time-scale of radial decay in exponentially thinning slender 
filaments of dilute polymer solutions. This relaxation time has been shown in simulations 
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to be proportional to the ratio Wi cs /Wi sc [24], while Clasen et al. [16] have experimentally 
observed that it is strongly enhanced with concentration even for solutions with c/c a st « 1. 
Firm confirmation may be obtained by more systematic rheological measurements [12] and 
by visualizing stained molecules in extensional flows [10]. The predictions may further be 
tested in simulations of multiple hydrodynamically interacting chains in elongational flows 
with appropriate periodic boundary conditions [17]. Such studies will also throw light on 
dynamical hydrodynamic and excluded volume screening in the dilute-to-concentrated so- 
lution crossover in polymer solutions. This understanding is critical to many applications 
such as turbulent drag reduction, ink-jet printing etc. where small quantities of polymeric 
additives are employed to induce large elastic stresses in the fluid in flows with significant 
extensional components. 

This work was supported by a CPU-time grant on the National Computational Infras- 
tructure at the Australian National University, Canberra. 
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I. CONFORMATION- AND CONCENTRATION-DEPENDENT FRICTION CO- 
EFFICIENT 



The friction coefficient C of partially unravelled and partially aligned chains at any concen- 
tration is anticipated to lie between the value Cc for isotropic coils at the same concentration 
at equilibrium, and an estimate Cr derived assuming that chains are slender rod-like objects 
of average length £ and d, all aligned in the principal stretching direction in a flow. A simple 
linear "mixing rule" is used for C, interpolating between Cc and Cr as follows: 

( z ~( z L-£ Cz^-V U 
where Cz is the limiting ( "Zimm" ) value for isolated coils at equilibrium, and L is the polymer 
contour length. 

Three distinct concentration regimes can be identified for modeling Cr/Cz depending on 
the instantaneous pervaded volume fraction of stretched molecules: 

where c is the molecular number density, c* is the conventional definition of the critical 
overlap concentration, which is estimated here as c* = d 3 for isotropic coils of diameter 
d = i?o/V3 at equilibrium; R is the root-mean-squared end-to-end distance at equilibrium. 



A. Cr/Cz in dilute regime, <j> < 10 



For non-overlapping chains, the mean transverse separation h is such that c£h 2 = 1, and 
hence from Eqn. (2), 

The following approximation is based on Batchelor's [1] results for suspensions of slender 
rods: 

K {L\ (4) 
Cz K+lnF \d a ) ' y ' 

where the argument of the logarithmic correction in the denominator is 

f £/d £/d 



H , if e/d > 1 + £/h 
1 . if £/d < 1 + £/h 



(5) 



2 



With this, Cr Cz as i -> rf , and for InF » K, 



Cr g I 

Cz ]n(£/d)-ln(l + e/h) d ' 



(6) 



The conditional specification of F in Eqn. (5) avoids singular behaviour as the aspect ratio 
£/d unity in the limit of infinite dilution, <fi -*■ 0; the choice of F = 1 is only needed for 
£/o? w 1 which typically occurs near equilibrium, where the mixing rule in Eqn. (1) ensures 
that Cr is of little consequence. The constant K = 0.15 is chosen to qualitatively match the 
hysteresis in steady-state extensional viscosity predicted by Brownian dynamics simulations 
in the infinite dilution limit (see Fig. 2 below). The approximation in Eqn. (6) is expected to 
be valid when h/d » 1 or yjcj) « 1 [1]. Therefore, Eqns. (4) and (5) are used when a/0 < 0.01. 

B. Cr/Cz in non-dilute regime, <p>l 



FIG. 1 Schematic of length-scales in a solution of partially unravelled, overlapping polymer chains 
aligned along the principal stretching direction in a uniaxial extensional flow 

Each molecule is a chain of N d = £/d "beads". Within each bead are correlation blobs 
[2; 3] of size £, with a total of N h blobs along the chain. Intramolecular hydrodynamic 
interactions are restricted to segments within a blob, but are screened at larger length 



d 
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scales. The hydrodynamic screening length £ h is that at which the intramolecular segmental 
density is equal to the average segmental density of the solution as a whole, which in turn 
implies that blobs of all molecules span the solution volume or 



cN h e h = l. (7) 

Assuming that an equilibrium-like structure survives at length scales smaller than the size 
of the transverse fluctuations d when Wi < 0(1), for solutions near the 0-state, 

f N K /N h ijd 

d* ~ N K /N d N h ' 1 ) 

where N K = L 2 /Rq is the number of Kuhn segments. The rod-like mean friction coefficient 
of the partially-stretched chain chain is calculated as the Rouse drag of iVh blobs, each with 
a Zimm friction Q = (£h/^o) Cz, 

S>Z "0 "0 

C Cr/Cz in transition regime, 0.01 < < 1 

Since screening sets in at incipient overlap, from Eqn. (9) above 

C d f 

= 7V d -y- = , at = 1. (10) 



Cc, z d do 

For a given c/c*, and given £ and d, if as calculated by Eqn. (2) is in the transition regime, 
Cr is calculated by linearly interpolating in between the value calculated above at = 1, 
and that calculated using Eqn. (4) with = 0.01; that is, 



Cz 



K 1-0 0-0.01 

+ 



K + In F|^ =0 01 1 - 0.01 1-0.01 
D. Concentration dependence of Cc/Cz 



The variation of Cc of equilibrium coils due to intermolecular hydrodynamic interactions 
in the dilute regime as c increases towards c* is not well studied. Beyond c*, the effects 
of hydrodynamical screening and crossover in dynamical properties in the linear response 
regime (e.g diffusivity) have been well documented [3]. These observations suggest that the 



variation in ( c with concentration may be approximated as follows: 

1 , if c < c* , 
Hoc*, 



Cz 



c 
c* 



(12) 



the second result following from applying blob arguments for equilibrium coils. 

As noted in the main body of the paper, screening is considered to be complete in either 
the stretched or coiled states when ( c or Cc reach the Rouse friction coefficient, estimated as 



Cz d Q 



(13) 



where b K = Rq/L is the Kuhn length. 



II. BROWNIAN DYNAMICS SIMULATIONS OF ISOLATED FLEXIBLE POLY- 
MERS 

Results of single-molecule Brownian dynamics simulations presented in this study are 
from Ref. [4]. The molecular model and its parameters, and the simulation algorithm have 
been reported in Refs. [5; 6]. Briefly, polymer molecules of N K Kuhn segments are mod- 
eled as chains of iV < N K beads connected by Finitely-Extensible Nonlinear Elastic (FENE) 
springs, and intramolecular hydrodynamic interactions are incorporated through Rotne- 
Prager-Yamakawa tensors. Excluded-volume interactions are neglected for chains at the 
9-state. Independent trajectories of chains are generated by integrating a stochastic differ- 
ential equations for bead positions that describe their motion under the combined action of 
hydrodynamic drag forces exerted by a homogeneous flow of the surrounding solvent and its 
thermal fluctuations, and intramolecular forces arising from free-energy interactions. The 
gyration tensor 

1 N 

G^EM), (14) 

where R u is the position vector of the v-th bead relative to the center-of-mass of a chain of 
N beads, and the angular brackets represent averaging over a large ensemble of independent 
stochastic realizations of chain configurations. The volume pervaded by a chain V = Vdet G. 
In flow, the polymer contribution to the extra-stress T p is given by the Kramers' equation 
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FIG. 2 Coil-stretch hysteresis in steady-state extensional viscosity: symbols are results of Brown- 
ian dynamics simulations with isolated 26-bead chains, obtained with initial ensembles consisting 
of equilibrium coils (A), and consisting of initially stretched chains (▼); continuous curves are 
predictions of the conformation-tensor model at infinite dilution. 



[7], and the extensional viscosity is denned as 



% = 



(15) 



The evaluation of the longest relaxation time for denning the Weissenberg number Wi is 
described in Ref. [6]. Predictions of coil-stretch hysteresis are obtained by running to 
steady-state two separate ensembles with different initial distributions of chain configurations 
at each strain-rate e: results for coiled states are obtained with an initial ensemble at 
equilibrium, while stretched state simulations start with chains nearly fully extended and 
all oriented in the ^-direction. Figure 2 shows hysteresis in fj p obtained with simulations, 
along with predictions of the conformation-tensor model in the limit of infinite dilution, with 
K = 0.15. 
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